function [V, D] = symeig(A)
    % 使矩阵对角化。
    % D = symeig(A); 返回与A相似的对角阵。
    % [V, D] = symeig(A); A = V*D*V^-1，其中D是对角阵。
    r = A.diagonalize();
    if nargout == 1
        V = r{1};
    elseif nargout == 2
        V = r{0};
        D = r{1};
    end
end